General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



JPL PUBLICATION 82-88 





An Investigation of Impurity 
Centers in Semiconductors 
of Variable Composition 

Part I: General Theory and Some Applications 


O.H. von Roos 



fNASfl-CE- I 6965S) AN IN VES 3IGATI0N OF 
IMPORIIY CENTERS IN SEMICONDUCTORS CF 
VARIABIE COMPOSITION. PARI I: GIKERAL 

THEORY AND SOME APPLICATIONS <Jet Propulsion 
lab.) 54 p HC A04/HF AO 1 CSCL 20L G3/76 

October 15, 1082 ' , 


N83- 15 148 


Unci as 
02250 


I Prepared for ^ . . • 

U.S. Air Force Wright Aeronautical Laboratories - 

' Wrrght-Patterson Air Force Base, Ohio 

Through an Agreenient with . . _ , \ 

National Aeronautics and Space Administration 

by ^ - ■ ■ ; ^ 

Jet propulsion Laboratory 
California Institute of Technology 
Pasadena, California' 


vim lun^i K^AiioN IV m 


An Investigation of Impurity 
Centers in Semiconductors 
of Variable Composition 

Part 1; Goneml Theory and Some Applications 
0,H. von Roos 


Odoboi U1. unv,' 


Piopivovi Km 

U.S Air Foigg WilQht AGitmmitiool Ubomtorios 
WilQhl*P(\l(Qr90n Air Foigg Bgsq, Ohio 

I hiouiih m Abtoomonl with 

Noltoiuii Aoronmilioa niui Spaoo Aclmintatrotloo 

by 

Jol Propulalon Uboratory 
Oolitomui liuUilulG ot Toehnology 
Pt\«Kdono. OnlilomiG 



Prepared by the Jet Propulsion Laboratory, California Institute of 
Technology, for the U.S. Air Force Wright Aeronautical Laboratories, 
Wright-Patterson Air Force Base, Ohio, through an agreement with the 
National Aeronautics and Space Administration. (NASA Task RE 182, 
Amend. 106, Air Force MIPR No. FY1455-82-N0607) . 



Abstract 


A theory of deep point defects imbedded in otherwise perfect semiconductor 
crystals is developed with the aid of pseudopotentials. The dominant short- 
range forces engendered by the impurity will be sufficiently weakened in all 
cases where the cancellation theorem of the pseudopotential formalism is oper- 
ative. Thus, effective-mass-like equations exhibiting local effective potentials 
derived from nonlocal pseudopotentials are shown to be valid for a large class 
of defects. A two-band secular determinant for the energy eigen-values of deep 
defects will also be derived from the set of integral equations which corresponds 
tj the set of differential equations of the effective-mass type. Subsequently, 
the theory in its simplest forl'^, will be applied to the system AlxCa^-x^s ; Se . 

It is shown that the one-electron donor level of Se within the forbidden gap 
of AlxGa]^_xAs as a function of the AlAs mole fraction x reaches its maximum of 
about 300 meV (as measured from the conduction band edge) at the cross-over 
from the direct to the indirect band-gap at x = 0.44 in agreement with recent 
experiments . 
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1 . 


Ivttroductlon. 0i\o of the outsCunding problems in modern day solid state 


physics continues to be the determination of the properties of "deep" 
impurities in otherwise perfect semiconductor crystals. By "deep" point 
impurities we mean either substitutional or interstitial foreign atoms 
embedded in the host lattice or just the absence of a lattice atom, a 
vacancy, with the proviso that the defect is capable of binding either one 
or more electrons or one or more holes in various energy states whicVi 
are lying deep inside the bond-gap. This distivigulshes deep impurities 
from donors and acceptors, also kno^«i as shallow impurities, their oxiergy 
levels being situated very close to a band edge. Whereas for some time now, 
the theoretical understanding of shallow levels has reached a quite satisfac- 
tory state based on the hydrogen mode, the diff lenities encountered in 
the study of deep impurities largely stem from the existence of a relatively 
strong and rather short-range potential which valence electrons experience 
in the vicinity of the impurity. This is particularly true for isoelec- 
Uronic impurities since there exists no Coulomb tail in this case, the 
potential being entirely confined to the central cell. The major difficulty 
presents itself in the fact that the extremely short-range nature of the 
interaction potential with its covicomitant abundance of high Fourier com- 
ponents in wave vector or crystal momentum space prevents the applicability 
of effective-mass theory used so successfully in dealing with shallow levels 
(donors and acceptors). Much progress has been achieved however during 
the last decade concerning the deep-impurity problem. Here is not the 
place to review past achievements in the theoretical understanding of deep 
impurities nor do we want to discuss in any detail the Importance of deep 
impurities as recombination centers governing the lifetime of electron- 
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hole pairs, a quantity so important for semiconductor device technology. 
Suffice it to cite recent review articles on the subject of both the 
theoretical understanding, or lack thereof, and the importance for device 
technology. The reader will also find a plethora of references on the 
subject in the literature cited. 

As we have mentioned earlier, it is the strong, short-range potential 
engendered by an impurity which constitutes one of the main difficulties 
in solving successfully the impurity problem. This situation may be 
ameliorated by the use of a pseudopotential formalism, and the aim of this 
report is to give an account of a theory of deep impurities encompassing 
such a formalism. To be sure, there have appeared in the literature a 
number of accounts on the application of pseudopotentials to the impurity 
problem. 3)22)34) The reason for the idea to employ pseudopotentials in 
the kind of problems we are discussing here lies in their very nature,^) 
Pseudopotentials and pseudo wave functions define a mathematical transfor- 
mation of the original Schrodinger equation which, while leaving Invariant 
the eigen-values of the energy, modify the wave function and the potential 
in a certain desirable manner. By orthogonalizing wave functions belonging 
to higher quantum states (valence bands) to the tightly bound inner electron 
shells (core states), this transformation achieves a weakening of the 
original potential in the inner core, precisely where the original 
potential, as seen by an outer electron, is strongest. This so-called 
cancellation, theorem has, as a consequence, the agreeable feature to 
concomitantly weaken the high Fourier components in a crystal momentum 
expansion of the perturbed impurity wave function, thus making it possible 
to use effective-mass theory in cases where hitherto it was thought to 
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be inapplicable. The ; -.oneering work of Pantelides3)22) bears ample witness 
to the truth of this conjectvire. 

In this report, we again tackle the problem of deep impurities with the aid 
of the pseudopotential formalism. In contradistinction to previous work3)22)^ 
we stay in the Wannier representation throughout. The advantages are twofold. 
Firstly, it is easier to ascertain the validity of -■ always inevitable - 
approximations made during the course of the derivation, and secondly, it will 
turn out that the equations governing the envelope functions peculiar to the 
impurity problem become local partial differential equations in contrast to 
the nonlocal equations derived previously .3) 

Section 2 of this report concerns Itself with pseudo-Wannier functions of 
the perfect crystal which form the basis for an expansion of the impurity wave 
function. Although pseudo wave functions and pseudopotentials are not unique 
generally, we confine ourselves to the use of the Austin form of the pseudo- 
potential^), in which case there exists a unique one-to-one correspondence 
between ordinary (true) wave functions and pseudo wave functions. The pseudo- 
Wannier functions will turn out be less localized than the ordinary Wannier 
functions, and to what degree this might be detrimental will be discussed. 

In section 3, effective-mass type equations for the motion of electrons in 
the presence of a point impurity will be derived. The role of screening of 
the bare electron-ion potential mediated by electron-electron interactions 
will also be discussed. The main part of section 4 deals with the derivation 
of a secular determinant for the energy eigen-values from the set of integral 
equations, corresponding to the differential equations found in section 3, 
by a method first Introduced by Bassani et al.29) Because of the use of 
pseudo-potentials , the number of bands to be considered in setting up the 
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secular determinant may be limited to but a few. In section 5 finally, the 
theory developed in section 4 will be applied to the system A& Ga. As:Se. 

Since selenium is a donor substituting for arsenic, an occupied level fairly 
close to the conduction band will be present so that the theory in its simplest 
form, taking into account only the condviction band with its various equivalent 
and non-equivalent minima without intervalley mixing, seems to be adequate to 
compute the donor energy level as a function of x, the mole fraction of AlAs 
in the compound under consideration. It is found, in agreement with experiment, 
that the magnitude of the energy level exhibits a maximum of '^J300 meV at 
X = 0.44, i.e. at the cross-over from the direct to the indirect band edge. 

The calculations reported here, have been performed by means of a variational 
principle and constitute one of the few cases in which computational labor 
could be held to a minimum. 
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2. Preiimlnarles 

Before we embark on the general theory of imperfections, let us review 
briefly the pseudo-potential formalism as ii applies to a perfect crystal. 

Let the Hamiltonian of the perfect crystal be 

H = T + V , (1) 

where T signifies the kinetic energy and V the periodic potential energy with 
the property 


V(H + “ VCr) 


U) 


Here ^ a^ , a Bravais lattice vector, is defined by the arbitrary 

j 

C^t') 

integers nj' and the primitive vectors a^ of the unit cell. For crystals con- 
sisting of a single kind of atom, for instance silicon, we have 


V 



v(r - R ) 
~a 


(3) 


where the sum runs over all sites within the periodicity volume V, and v(r) 
signifies the potential energy of an individual ion located at the site R^. For 
crystals containing more than one kind of ion within the unit cell, for instance 
GaAs, we may write 


V = 




(4) 
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In this case, the additional sum runs over the various ions j located at positions 
b, within the unit cell defined by the. lattice vector R . With the definition 

J •'U 

v(r - R ) o \ V. (r - R^ - b. 

eqs. (3) and (4) become identical, and no mention need be made of the complica- 
tion arising from having more than one species of ions within the unit cell, at 
least as far as the development of the general theory is concerned. 

Schrodinger' s equation now reads 

(T + . C6) 

The Bloch functions, 



♦„k " ■ -%k<f) 


(7) 


therefore possess the energy eigen-values E^Ck), with k the crystal wave vector 

and n the band index in customary fashion. The quantity u^j^ is periodic in R^ and 

therefore satisfies an equation identical to eq. (2). It is to be stressed that 

2 ) 

wo imply the Hartree-Fock one-electron approximation when writing eq. (6) . The 

o) 

potential V is therefore, by necessity non-local. “ But these complications do 
not concern us here, since we assume that the unperturbed crystal eigen-values, 
E^^(k), are given, known functions of their parameters. 

For reasons which will become apparent later, we will transcribe eq. (6) into 
the pseudo-potential formalism in the Wannier representation (simiJar calculations 
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3) 

have been performed by Pant'- tldes ). The reason for "pseudizing" che potential 
(3) are amply described in ref. 4. Here we can only give a very brief outline. 
Suppose we separate the eigen-functions, into those belonging to the core 

•y 

5) 

states, and those belonging to the valence and conduction bands, * . The 
inner core electrons are usually strongly localized at the atomic sites whereas 
the outer valence electrons form broad bands and are therefore delocalized to a 
greater or lesser extent. To each wave function, satisfying eq, (6), we now 
associate uniquely a pseudo wave function satisfying 


(T + V)(}) 


vk 


-V 




v!s 


^vk 


( 8 ) 


The eigen-value, E^(k), of eq. (8) is identical with the corresponding one for 
corresponding to eq. (6). Furthermore, the connection betweeti and is given by 





(9) 


where the coefficients, 3 (k), satisfy 

C‘.V 


- 




6 (k) 

cv 




',>3 ' 

k c 


(k) = . (10) 


vk 


We have used the familiar Dirac notation throughout. As can be seen from 

6 ') 

eq. (8), we have used the Austin form for the pseudo potential. ' From eq. (9), 
it follows that 


“ "*cki*vk" 


( 11 ) 


7 




because and are orthogonal, Sometimes the argument is turned around by 
writing 


^vk ” *^vk ^ *^^ck^ ^vk^’^ck 

c 


( 12 ) 


which shows that with any choice of the wave function is made orthogonal 
to the core states by construction. This is important for practical calcula- 
tions where convergence of a secular determinant toward a vrlence band energy 
eigen-value rather than toward a core state is desired or where a variational 
principle is desired to "home in" on a valence band energy rather than on a core 
state energy, the latter being several keV below the valence-band energies. As 
an added bonus, the convergence of the secular determinant in a plane-wave expan- 
sion, v^^'-g simple plane waves to represent the wave functions becomes much 

more rapid when using the device of eq. (12). Known by the name of the orthogonal- 
ized plane-wave method (OPW method) in this case, it had been proposed by Herring 
long before pseudopotentials became fashionable. ■' ' 

9) 

As seen from eq. (8), the pseudopotential is nonlocal , but we will see 

later that, nevertheless, the impurity problem formulated here within the 

pseudopotential approach will give rise to a toaat effective potential. The true 

advantage of the Austin pseudopotential employed here consists of the validity of 

the cancellation theorem. This theorem states that the pseudopotential of 

11 ) 

eq. (8), which we may rewrite concisely as 


V 


A 


V 





(13) 
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Ih I’onKldin'ahly wonlcor Xn the eore rogloiiB than the original potential V wlileh 
la rather atrong there. “ Thla la of paramount Importaneo for the impurity 
problem, aa clearly realiKeil aome time ago. But before going Into this 
matter, let us trimserlbe eqs. (8) and (10) Into the Wannler representation. 

The Wannler functloua are defined by 


rl/2 


-Ik* R 


w.,(n) « N ^ e ' 


uk 


(l^'ta) 


wllli the Inverse transformation being 


• 1 /: 




- N • e ' V(a) 


nk 


(I4b) 


Thi' normal Ir.aL ion adopted lu're Is the following; • v , |m « vS' ,iS, i, 

nk'Mi k nn kjk 


/. f 


Ik-R 


\ ’ 


ik*R 


■ a 


E 


^ N 7 'V|k ’ N Is tlie number of unit 


ci'lls In the quantisation volume V, is a reciprocal lattice vector and 

1 0) 

sums ovi'i* k extend over the first Brillouln xone. Also, for a poten- 

tial as given by oc], (4), the expression (14b), for instance, becomes 


'^nk 


N-iC.: Vo'"’"" 


V Wj^'^^(r “ - h j). For simplicity of notation, we 


write w (.t) " w (r “ R ) 
n n a 




' » c 


9 



ORlGlNAt PAGS 
OF POOR QUAV.1TY 

Analogoufj to imih. UA), tUo Wanuloi’ pmnulofuiu’t lonts wav bo lutrodnood via 

~lk‘R 
k 

ami 

lk*U 

Z-j"' 

a 

Nolo that tUo band tndox n of oii». (lA) onoompanaon all poaHibnitioa fooi'o and 
valonoo alatoa'I whoi'oaa tho band indox v in oqa. U?>) Xu I’ouorved to valonco 
H La I on only, 

InunHlnolnp, now ibo Ikniolor tranuloian oi tbo ooolfto lontu, , doflnod In 

oq. (101 a. a 


W^(a) 


N 


■I/: 


/ j 


with tho invoTvai' 


d (vil c! N 
ov 




“Ik’R 

,(k)o 


ik’H 


E A »V 4\ 


(Kia) 


(Ibbl 


wo obtain from oqu. (0) avul (lb) tho oonnootlon botwoon tho puoudo-Wannior fnuoiionu 

1 ^ 

W and tbo rogular Wa;ml or fvnu'tlonu w , namolv tbo rolatlonuhlp 


W (a) 

V 


w (a) 
V 


>T~A 

/ 1 


0 , \ 


(a 

0 V 


Ow. (\) 


(in 
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llu' WauuU'V luuv‘1 hM\a, «uv U'chUiuhI nhoul Uu'iv roniunH ive alomU' 

1 }") 

at Ion, wo noUoo that I ho i>ao\ulo K^mulov tuuot loua, uro UvalL^oU alnnu 

i\ll niton nn Iho muu ovor t tu oq, lovoatn, lUnv’ovor, thta lf4 aooomt'l tnlunl 

with vnn*tv\j\ auoooaa Uo|\oiuHu}\ on tho watwttmlo v't tho oooUlotoula ‘^n 

n luwot ton ol \ Iv^r l tnovt a, \w ovvUm' to trattnovlho oqn. tiO uml into tho 

W.mntov voproaontiU ton, wo wlU omplov tlto lot lowing vlovtoo, Hootmno tho j'oton 
tlitl V !n in'otodio toq, tOl, wo hnvo in qonooal thvit 


nl\ 


,lv| 


nlv 


vN, 


k’jlv'^u'K 


V| 


ulv 


vuo 


Thin in tho vonnon whv, In oq, ttil, tho nniu ovor tho ooro at at on vnttv oxtonUn ovor 
tho tnunl tnvlox o,Kooi>lnH K 1 Ixo^l. \vV now vola\ thin voqnt routont . V'qnnt Ivma tti) 
amt ttion f.o ovor Into 


cr t v'tv" , 7 1 ''V , - i- , 

/ i vlv ok ' * vk Y vk 
v‘k ’ 




attvt 


'*\*k 


w<sesc3s« 

'vk * / J ^''ok 

Ok’ 


U'O't 


Fqttatlon tit''' iuH'omon analoy\vntaiv 




OmIv 


"t 




0 V 


tk,k 
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Equations (19) and (21) are, of course, equivalent to eqs. (8) to (10) provided 

19) 

that eq. (18) is used. Now the transformation into the Wannier representation 

with the aid of eqs. (14) and (15) becomes easy. But first let us define the 

1 7 18 

Fourier transforms of the band energies, E^(k) » 


coiicomitantly with 


^ (a) = N 
n 


-1 




-ik-R 


(k)e 


__ lk*R 

E, (k) = > <? (<x)e ■ 

1 1 X i 1 L 


(22a) 


(22b) 


Introducing the Fourier coefficient, 


, v-^ -ik-R +ik'-R 

2-. ' ■" ' ■' - 

k,k' 


(23) 


obtain from eq. (21), after some algebra, 


E 


^^(^s - Y)i:^^^(a,.S) - ^^(u - cn i\,^(S,l) 


- / ^ <w^,(\) jvlw^^.(vS) ' = 


c',6 


(l) |vlw^(u)r' 
(24) 
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If WG now restrict ourselves to the symmetry of eq. (18), becomes only a 

function of the difference (a-8). In this case, eq. (23) goes over into eq. (16a), 
and equation (24) becomes 





■) 


3 (a 

cv' 


Y - 6) 


-E 

c'6' 


<w^(y) 


w ,(6)>3 , 

I r ^ , 


(a - 6) 


<w^(y) | viw^(a)> • 

(25) 


It is eq. (25) x\rhich is the counterpart of eq. (10) in the Wannier representation. 
A direct transformation of eq. (10) with the aid of eqs. (14) would have met with 
difficulties. This is why we first relaxed the symmetry condition of eq. (18) aad 
imposed it after the transformation luad been performed. 

The transformation of eq. (19) into the Wannier representation follows along 
the same lines and is given by 


(T + 


V)W^(a) 


-E 

c,y 


Wp(Y) 


'■"c 


V 


W ((x): 
v' 


=E 


S' (y - 

V 


e)W^(y) 


(26) 


Before x^fe go to the heart of the matter, the impurity problem, let us 

Irivestigate the meaning of eq. (17). The pseudo-Wannier function, eq. (17), 

evidently is localized about all sites, as already mentioned. The coefficients, 

3 , wlilch determine the localization about centers other than the main site a, obey 

cv’ 

eq. (25). Now, in order to deal with impurities, particularly deep centers with 
their strong localized potentials, one wishes to employ wave functions which are 
themselves localized so that a perturbation scheme of one type or another may 
converge rapidly (it is hoped). At the same time "pseudizing" the potential 
weakens its central core, so that high Fourier components in k-space become small, 
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and approximations patterned after the effective-mass approximation may become 
acceptable even for deep centers (bound energy states deep in the forbidden gap 
between valence and conduction bands). But at the same time, the pseudo wave func- 
tion of eq. (17) becomes somewhat delocalized. How much delocalization takes place 
is governed by the coefficients 3^^. Because of the strong localization of the 
core Wannier functions w^(y), we have 


<w (y) Ivjw ^(a)> « V 
c c 

to a very good approximation. Putting a - y = 

<w (y) |v| w (a)> = V 

c ' ' V 


.5 

cc ya 


20 ) , 

t; and 


cv 


(?) 


(27) 


(23) 


as well as 




(29) 


equation (25) goes over into 

Y c.' 

Although it is difficult to solve eq. (30) exactly, it is quite easy to See, at 
least qualitatively, that the dependence of 3^^(?) on Cj the distance between any 
two lattice sites, is governed by 

equal to zero, the only solutions for the coefficients 3 are since the 

determinant of the coefficients of the ensuing linear equation for the quantities 

21 ) 

3 (?) does not vanish in general. Introducing the potential matrix element 

cv 

V (c), we see that 3 U) will become proportional to V (?) and in higher order, 
cv cv uv 
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will become a function of V V etc* But from the definition V (^) of eq, (28), 
we note that, due to the strong localization of the core state w^(y), the matrix 
element will exponentially decay with increasing G* Therefore, the 

delocalization inherent in eq. (17) is rather weak and may be ignored in many 
applications . 

3 . The Point I inn urity Problem 

Suppose now that the perfect crystal contemplated in the previous section 

contains an isolated impurity. This impurity may be a substitutional or 

Interstitial atom or it may be a vacancy. Whatever it is, it will give rise to 

an additional potential energy which we denote by U. Rather detailed analyses 

of the structure of U have been presented by several authors, particularly for 

21^3)22)23) fairly straightforward to generalize this analysis to other 

semiconductors. But we will treat U for the. time being as a given function of 

all its relevant parameters, in particular its spatial coordinates. In order that 

no confusion arises, let us introduce the following notation. All unperturbed 

quantities (U ~ 0) will be annotated with a superscript zero. Therefore the 

0 

periodic potential V of the previous section now becomes V , the Wannler function 

w (a) becomes w^(a), the new quantities remain unprimecl, etc. Also, since the 
n n 

implantation of a foreign atom (or removal of a lattice atom) will alter the 
equilibrium positions of at least the adjacent atoms from to R^ say, we 
denote this fact by a primed Greek letter. A core function w^(a) of the 
unperturbed lattice now becomes w^(a'), for instance. We assume that the imper- 
fection is located at a = 0 (R^ = 0) or, for an interstitial atom, in the immediate 

neighborhood of « = 0. Incidentally, the collection of quantum numbers c for the 
core states differs in general for foreign, atoms from those of the host atom 

except of course for isc/coric atoms. This fact must be duly taken into account, 

24) 

but we will not mention it explicitly in the ensuing formalism. 
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With = T + V^, Schrodinger' s equation in the "pseudized” version for the 
system-perfect lattice-plus-impurity - reads now, in analogy to eq. (26), 

(H° + U)ijj - w^(y")<w^(t") |'V° + U|ilj> = , (31) 

c,y" 

where the sum runs over all core states and all sites y'. We propose to solve 
eq. (31) with the following ansatz: 


'i' = ^ F^(a)W®(a) 

v,a 


(32) 


We note that the sum extends only over valence (and conduction) bands, and the 
positions a are taken as those of the unperturbed crystal. This latter fact 
does not constitute an approximation, since the envelope functions will depend 
parametrically on the a". The pseudo-Wannier functions W^ do not form a complete 
set unlike the ordinary Wannier functions, w^. Had we expanded tjj into the com- 
plete set w^, we would have had to include core states n = c in the expansion 
even in the limit U = 0, since for U = 0, eq. (31) goes over into eq. (8) 

(or eq. (19)). In this case, any of the W^ themselves become solutions of eq. (31), 
and according to eq. (17), even then core states are present. In other words, there 
would be no hope to confine an expansion of into regular Wannier functions 

n 

tVi 

to include only a few valence bands, since even in 0 order, ip would already con- 
tain a sizable emiount of core levels. The ansatz (32) truly represents an approxi- 
mation in anticipation of .he need for only a few terms in the sum over the bands. 
From the orthonormality of the Wannier functions^^\ we have 


<w^(a) |w^^(B)> 
n n 


6 

nn 



(33) 
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and from eq. (17), it: follows that 




( 34 ) 


Multiplying eq. (31) by <w^(a) |, using eqs. (22) and the wcll~kno\m relation^®^ 


<w^(«) !h^|w^.(3)^ = - «) 


(35) 


we obtain a set of equations for the envelope functions F^(a): 




'< - W^, (D I v° + U 1 > ( tV) , ( p) == EF^ (a) . 


(36) 


Since the core states are well localized, we make, the, further approximations 


:w^(a)jw (r')'’^ .<w^(oi)jw (a')'' 

V ‘ c n 't V ' c 


(37a) 


and 

M^r^,(Div° + UjW°.(P)'- ^ '%(^^')1V° + ulW^.(P)'' . (37b) 

The meaning of these approximations is given in the following. The equilibrium 
positions R of the host lattice are slightly displaced to R' = R + ^R due to the 
perturbation U. Owing to the strong localization of the core states about their 
respective sites and also to a lesser degree the localization of the Wannier func- 
tions for valence bands around their own sites, the matrix elements (37) are 
expected to be largest when the sites almost coincide. We presuppose that 

/R 1, so that there is no intermixing of sites, 
a >x 
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Next we consider the matrix element, 

= <w°(«)|u|w°.(0)> . (38) 

Away from the impurity (we remember that the impurity is located at a = 0 
(R^ =0)) the potential U is smooth, so we may approximate'*' ** we have 

In the vicinity of the impurity where U varies rather strongly, this approximation 
is not allowed and the matrix element must remain as it is. Equation (39) tells us 
that the second sum over 8 on the LHS of eq. (36) reduces to a single term if a is 
sufficiently far away from 0 (the position of the impurity). 

Finally, with the well-known expression 

^ ^^(3 - a)F^(3) = E^(-i7^)F^(a) , (40) 

3 

17) ] 8) 

(which follows from eqs. (22) together with ' ’ , 

R .V 

e ^ ^F(y) = F(y + a) , (41) 

where 7^ signifies the gradient with respect to a (or R^)) we obtain the follow- 
ing set of equations for the envelope functions from eq. (36): 


E^(-iV^)F^(ot) +Y\ %(a)| TJ 


V 


w (a^)><w (a'') I V + U 
c c ‘ 


W .(a)>F .(a) = EF (a) 

V V V 


(42) 


18 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Before eq. (42) can be put to use, the potential U must be specified. But first, 
some general properties of eqs. (42) which are independent of the precise form 
of U may be pointed out. If a (or R^) is considered a continuous variable, as 
is customary eqs. (42) constitute a set of lochil partial differential 
equations, although the underlying theory is essentially nonlocal. For an Iso- 
coric impurity, wo hove 

w^(a') « w^(a) (43) 

for all a, even in the neighborhood of the impurity coll, and eqs. (42) go over 
into 


v' 

owing to the orthogonality of the valence functions w^ with the core states 
Furthermore, far enough from the impurity wheiT a becomes equal to a', the core 
states w^ also become equal to w^, and again eqs. (44) become valid. This is true 
whether the impurity core states are isocoric with the host core states or not. 
Therefore, the pseudopotential of the interaction matrix elemcirt in eq. (42) 
becomes effective only within the central cell of the impurity center and will 
effectively considerably weaken the ii\teraction potential U precisely there 
where it is strongest. This is another maviifestation of the cancellation 
theorem. Similar findings have also been reported by Pantelides. 

The Interaction potential energy U consists of a number of terms. Signify- 
ing the interaction of the valence electrons with the rest of the system and with 
each other, it consists of the bare interaction with the nuclei and all the core 




19 


original B 

OF POOR QUAUT'r 

. 25) 

states U^, together with the oolf-aonQistcnt electron’-electron interaction. In 
the one-electron approximation in which an electron satisfies the Schrodinger 
equation, the analog of the "pseudized" version of eq. (31) is 

(H° + ■= Etj; , U = (45) 

in which the potential U is given by the sum of the bare (or external) potential 
and a polarization or screening potential U^. This latter potential can be 
thought of as being due to the readjustment of the electron configuration which 
existed prior to the application of the external potential into a new configura- 
tion. In the simplest case, the Hartree approximation, U^, is given by 

= e^y' d\^|r - r1‘^(lt|i(rOl^ - k^(r')l^j , (46) 

where i[i is the solution of eq. (45) and signifies the solution of cq. (45) with 
= 0. The solution of the non-linear system of eqs. (45), (46) constitutes what 
is generally known as the self-consistent solution and the potential U=U^+U^ as the 
self-consistent v^’tential. The derivation just sketched, when followed through in 

0(3 ■) 

detail, taking due account of the occupancy of all quantum states*" , results for 
each Fourier component of the bare potential in a correspondingly screened 


potential 
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The Fourier transform of U^(r) Is defined by Ujj(g) ° V d^vU^(r)e’'^3'£ with 


the 


inverse U, (r) 
b ~ 


E“> 


(q)e 


iq* r 


where the sum over q extends over all of q 


spare, not just the 1 ®^ Brillouin zone. The choice of discrete values for g 
Is consistent with the normalization imposed by the use of a periodic quanti- 
zation volume V throughout this work. 

A particularly simple expression for the dielectric constant e is obtained if 

S tl 

a) eqs. (45) and (46) are solved in 1 order perturbation theory (linearization), 

27) 

and b) if Umklapp processes are neglected. It is given by 


e(q) = 1 - ^ 
Vq‘ 




kv, v' 




F (k + q) - F .(k) 


V 


E (k + q) - E >(k) 
V - 2 V ~ 


(48) 


In eq. (48), the quantities F^(k) are the occupation numbers (1 or 0) of the 
unperturbed quantum states v, k and the E^(k) are the usual valence and conduc- 
tion band energies. The neglect of Umklapp or central field corrections seems to 

be justified if the bands involved are sufficiently broads Explicit calcula- 
28 ) 

11008“ have shown that Umklapp contributions are indeed small in most cases of 
29) 

interest. “ Underlying expression (48) for the dielectric constant is yet 
another approximation t<? all those already made. It is the assumption that 
the polarizability of the core states may also be neglected. This is expressed 
by the fact that the summations in eq. (48) extend over valence bands only. But 
this neglect is justified since the core states are tightly bound and therefore 
much less susceptible to the influence from the outer electrons. To summarize, the 
dielectric constant given by eq, (48) has been obtained under the stated 
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approximations by solving the self-consistent eqs. (45) and (46), Equations (48), 

being based on the Hartree approximation, by necessity ignores exchange and corre- 

25) 30) 

lation contributions , but these contributions give rise to small effects 

which we ignore. 

One more remark is in order. The analysis for the screening of the bare 
potential outlined above is based on eq. (45) which is the analog of eq, (31) 
but does not contain the pseudopotential Included in eq. (31). The pseudo- 
part of the potential becomes significant only within the central cell at the 
locations of the lattice and impurity sites owing to the presence of the core 
states w^(a^). It is precisely in these regions near the ionic sites where 
Unklapp, as far as the dielectric constant is concerned, becomes important. 

But we have seen that Umklapp is negligible for the broad valence bands we are 
considering here. Therefore, as long as we are justified in using expression (48) 
for the dielectric constant, we do not have to screen the bare potential at 
all as far as its appearance inside the nonlocal pseudo-part of the total 
potential is concerned. We may therefore write for the total potential 


V . 
w 


(aO 


<w (a) 

V 


U 


E 


w (a')> <w (a"”) V 
c c 


+ U, 


I W^^ (a)> 

I Y 


(49) 


where 


U 



q 
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(50) 
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is the sereened bare potential with e given by eq. (A8). Equations (/i2) for the 

I 

envelope functions F^(ot) now become 


'■v'' 


■1V„)F^(«) 

v' 


V .(a')F„«(it) 

VV V 


KFy(a) 


(51) 


Equations (/i9) to (51) form the basis of our further investigations. Inasmuch 
as the dielectric screening, eq, (68), follows from I hu'ap response theory, the 

•V 

potential V must be considered to be weak, but this is true for most applications. 
The potential within the inner core is weak because of the cancellation theorem, 
and far away from the Impurity V, becomes Just 


1 im 
R“>»' 


V 


2 

t.(0)R VV 


(52) 


with the net atomic number, 2, and the static dielectric constant, » (0) . Then 

31 .) 

for small enough E, the total potc.itlal V is indeed weak. 

^ * Met hods o f S o 1 u _t i mi. 

There exist essentially two methods to .solve eqa. (51) for tlie envelope 

function F^(a) and the associated energy eigenvalues. Tlie first method consists 

of considering eqa. (51) ns differential eq.! .tions in the now continuous variable 

a (or R^ R) and expanding the. operator of the band energy E^(-iV^^) In a power 

series in ~iV^^ about an energy extremum ami breaking off this series usually 

after the quadratic term. This is called the effective mass approximation (EMA) 

and an excellent review of it is given by Stoneham. Complications arise \dien 

the bands are degenerate, which is usually true for the valence, bands Just below 

.3’’) 33) 

the conduction band in semiconductors of group IV and IXI-V compovmds. ‘ The 
FJ^iA is valid only for shallow states (donors or acceptors) when a ov\e-band 
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approximation suffices. However, it is also valid for many deep levels provided 
that the band gap is indirect, as shown in ref. 2 and 13. For semiconductors with 
a direct band gap, the one-band approximation is almost certainly incorrect for 
deep levels even though the effective interaction (49) has been considerably dimin- 
ished in the central impurity cell by the pseudopotential device. An attempt along 

34) 

these lines taking several bands into account has been made some time ago. Butin 

ft ^ 

this work, the Kleinman-Phillips version of the pseudopotential ' was used and the 
sum over the core states, the analog of our eq. (31), included also valence states, 
and no screening was taken into consideration. That theory is therefore dif- 
ficult to compare with eqs. (51). Before we leave a discussion of the EMA, we 
like to point out one simple approach which is applicable to eqs. (51) wlien they 
collapse into a single equation (one-band approximation) . We have seen in the 
last section, that the potential V^^^(R) becomes a slowly varying function of R 
(see for instance eq. (52)) for large distances from the impurity, whereas near 
the impurity the pseudopart of the potential avoids much of the variation 
of the original potential. For a single donor-like impurity, as an example, 
we m.ay therefore approximate the potential for the ground state (S-state) by 

I -e^/c(0)R for R > Rq 

(53) 

- ^ for R < Rq 

35) 

where ^ xs a suitably adjusted constant. The equation for the envelope func- 
tion F^(r) (the subscript c indicates that the level to be determined 
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>v 

»plU.s oil I ho ooiuluotion buml) , iuiBuwiuR a acalmr effocdive nuiaa i« foi* 
H I mp 1 1 city, b oc o mo 55 


h* ■' 

** 

»v ' 




K^.) l'V(U) 




K^. ia iho oiuM'gv ol Iho odpo ol Iho oomluolUm baiul and V(lO Us p,tvon by oq. ('ii), 
riu* bomularv oinulllUma, I.0,, p.ooa to koto at InlinUv and F^, sstayM finito at 
v'.oro lopolUor with tho ooiuluullv of iho loparllhmio dorlvatlvo ;u R R^, docor~ 
mlno lh»' oij’iHV'Valno lor F, oiuh' and aro known. Tho potontlal oonntt- 

Inlo.a a modo I potont lal, and an auoh ban boon InU’oduood, albeit in a somo.whaL 
vitt toronl oonioxi, lonp apo- Unnally, the oonatanta R and aro dotormlnod 

l oom oxporinunual dala, U wo ran adopt a "muff in~tin" model, thou the ohoioo of R 

I/) 

la fjvnnowhal arhll rarv and mav ho tlxod onoo and lor all. U\ft/ Is dotormhiod 

ns lap Iho oxpor imonU\ 1 I v knovm provmd state onorpy, all oxoitod states 

and iholr onorpios oan ho oompnlod uslnp iho sa?';t‘ oonstaut This lolUns’s from 

oq. 1>VU ond Is iiniti' a dlstlnotion from tho UvSnal monel potential for un)H’r-’ 

tnrht'd orvfstvils \vi\oro lor oaoh anpnlar momontum \',a now oonstaut lias to ho 

dotormlnod. 

For doop lovt'ls In d Ironmnind-pap somieonduotors, the ono-band FMA eonsti- 
lul os at host .1 ]'oor approximat Ion. Wo aro looking, therefore, for an alternative 
wav ol solving oi[s. tSll. Now, 11 tho onorpv oipou“Valuos of tho impurity are 
oxpootod to lie in mid”*band"-gap , then not only do the condnetion and valence 
hands oontrihnto, hut also oonsldorahlo parts in k-spaoo of tho hand onorglos 
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ivro involvtHl. To take full advantage of the situation, a Fourier trans 
fornualon of eqs. (‘il) is Indicated. Therefore, defining 




P (R ) 

V -"Ot 


■= N 


- 1/2 



(55a) 


£„0c) 



(55b) 


as well ns 


38) 


-ik-R 

” , (56) 

a 

we obtain from eqs. (51), 



5 - K^(R)] ,r^(10 V„„.(k - k-)f„.(lc-) . (57) 

v,'r 

The summations over k or k' in the above equations only extend over the 3®^ 
Rriiiouin sone. The integral equations (57) are equivalent to eqs. (51). 

Before going any further, we allow the quantization volume V to become infinite. 
Tn this case, we are dealing with continuous variables, both in k as well as in R 
space The transcription follows the usual rules: 




(58a) 
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where 


V = • <92 ’< 93) 


(58b) 


is the volume of a unit cell with the primitive lattice vectors a 
We also have V = Nv. The Kronecker symbols go over into their 
Dirac iS-f unctions 


1 j ^3 * 

respective 



V 


6(k - kO 



= V 6(R - R') 


Equation (56) now becomes 


(58c) 


V 

vv 


(10 



3 

d R V 


-ik* R 


vv 


,(R) 0 ~ 


(59) 


The integral equations (57) are solved readily if the kernel V(lc - k'*) may be 

written as a product V^(k)V 2 (k"). This can be achieved by using a complete set 

39) 

of normalized functions such that 



(60) 


(61a) 


(61b) 
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It is to be noted that ^ is not the complex conjugate to because the 
potential V may become negative or positive as a function of R depending on 
position. This is so because, for instance,, a totally attractive true potential 
(negative for all R) will acquire positive parts for small R due to "pseudizing. 
We also have assumed the interaction potential V^^^(R) to be centrally symmetric. 
From eqs. (57), (59), (60) and (61), it follows that 




.(k) 

VV 


v,m 


(2n)- 




( 62 ) 


Introducing the quantities 


im 

v^v|vv' 



Ck) 



(63a) 


and 


^m 

VV 



_m 
^ VV 


.(k)f^.(k) 


(63b) 


\3e. obtain the linear set of equations 


a' .y G*”i a" , 


(64) 


We now consider only two bands, assuming that the energy eigen-values within the 
band gap between the conduction band and the valence band are mainly determined 
by these two bands. We ignore the complications introduced by degeneracy of sev- 
eral valence bands. Designating the conduction band by v = 1 and the valence 
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bnr.-l by V = 2, we obtain four equations (they are still matrix equations in 
the space generated by the which may be written in matrix notation 


^’11 1 11 ^11 ^’ll|l2 ^12 

^21 “ ^’ 21,11 ^11 ^’21 1 12 ^^12 


^12 “ ^’12 1 21 ^21 **' ^’12 1 22 ^22 'j 

^22 “ ^ 22|21 ^21 ^ 22|22 ^22 j 


.( 65 ) 


The determinant of the coefficients of this linear system must vanish, and it is 
this condition which determines the possible energy values within the band gap. 
With the help of the defining equations (59) and (63a), we also have: 






^\jl/lk) 


In 


■E*;: 


Im 

]/. 


m 


inn 

^’jl/lk 


d-^kd^c' 




Ik 


(k') 


(2n) 


(H - 




k-)\ 


- k) 




In 


ij/lk 


(66) 


After some manipulations, eqs. (65) may now be reduced to (again in matrix 
no tat Ion) 


G. , I . , + . t , - l\ A, , + I _ A_ = 0 


11 11 11 12 1 1211 11 11 12 1 12 22 22 


and 


(67a) 


‘•’ 22121*^^2 ^ 2 lill ^11 (^ 22 i 22 ■’* ^ 22 l 21 ^^2 '*^ 2 l |22 -^^^22 ' 


(67b) 
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Here we have introduced the inverse matrices 33 and 33 ^ defined by: 


and 


^1 " ^ ” ^^12 1 12 


38 =. I - , 

2 21 21 


(68a) 

(68b) 


Eliminating 


from eq. 


(67b) gives 


^11 " *^21 111 



so that inserting this into eq. (67a), it is seen that the energy eigen- 
values of the impurities are determined from the condition 


(69) 



,-l 

’lljl2^-^l *^12 1 


1 + G_ t _ £55 . 1 1^ ^ ■’ 


2l|ll f^^2l|22 ^2 ^22 1 21 ^^22 I 22 


-d) 


^11 1 12 ‘^1 ^12122'" ^ * 


(70) 


Here the symbol DET signifies a determinant in the matrix space generated by 
the functions g introduced via eq. (60). Since the functions g (R) are largely 
arbitrary, a judicious choice may reduce the dimension of the (infinite) deter- 
minant of eq. (70) to a manageable, size. In principle then, it is possible to com- 
pute the energy eigen-values of a point impurity, be it a substitutional or inter- 
stitial impurity or even a vacancy, with the aid of eq. (70) provided that the 
underlying band structure of the host crystal and the potential energy of the 
defect are knoma and provided that a number of approximations made in the course 
of this analysis are justified. It is perhaps worthwhile in concluding this work 
to point out the various approximations which led to the basic equations (49) 
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through (51) and aliso to the determinantal eq. (70). The underlying, unperturbed 
crystal structure has been formulated li\ the one-electron approximatloii. This 
in itself does not constitute a limitation, since it is not necessarily predi- 
cated on the use of a single Slater determinant as it is in the. Hartree-Fock 
approximation. Correlations may well have been taken into account. We left 
this question largely open, assuming that the unperturbed crystal properties 
(its wave fmictions and energy bands) are suff Icieiitly well known in order to 
attack successfully the point impurity problem. The approximations made during 
the course of this work may be discussed in their order of appearance. First., we 
expanded the wave function i|j of the system - crystal + Impurity - in pseudowave 
functions i|)^ belonging only to valence and conduction bands, ignoring the core 
states. But the pseudo-functions contain an admixture of core states in a natural 
way. As far as only valence states are concerned, they form a complete, 
albeit non-orthogonal set since there is a one-to-one correspondence, barring 
accidental degeneracy, between and Furthermore, the pseudopotential 

weakens the strength of the true potential, as we have seen, so that high Fourier 
components and therefore admixtures of core states in the expansion (32) tend to be 
de-emphasized. Next, the approximations eqs. (37a) and (37b) diagonalize the cor- 
responding matrix elements. While this constitutes an excellent approximation 
as far as eq . (37a) is concerned, it is less so for eq. (37b). The culprit here 
appears to be the pseudo-Wann .er fuiictlon, W^, which is less localized than its 
counterpart, w^. But we have seen from eq. 30 and the discussion following it, 
that the delocalization is not severe, and so the approximation of eq. (37b) is still 
quite justified. Another approximation which has not really been used, however, 
at least as far as the general formulation culminating in eqs. (A9) to (51) is 
concerned, appears via eq. (39). It is self-evident, and it is certainly true 
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far away from the Impurity. The last covey of approximations to be discussed here 
briefly (since it has been done extensively in section 2) concerns itself with 
the screening of the bare potential or the Influence of many-body electron- 
electron interactions on the one-electron Hamiltonian. The most severe limita- 
tion in this case consisted in using linear response theory, the bare potential 
being screened by means of a dielectric constant as in eq. (47). But this implies 
that the potential U^(R) is to be considered weak iii the sense of first order 
perturbation theory. While this is justified in most cases of interest for 
large R, first-order perturbation theory will simply not do for small K. How- 
ever by means of "pseudizing" , the true potential had been weakened for small R 
to such an extent that first-order theory became viable throughout the wliole 
range of R. The neglect of Umklapp (central field corrections), core polariza- 
tions, exchange and correlation corrections to the dielectric constant and the 
neglect of the screening of the pseudo-part (the nonlocal part) of the pseudo- 
potential constitute additional approximations. They have been discussed in 
section 2 and their effects found to be small in many cases of interest. 

Finally, turning to the derivation of eq. (70), two additional approxima- 
tions have been introduced. Confining the expansion (eq. 32) to two bands only, 
constitutes one, and to assume the potential V to be centrally symmetric constitutes 
the other approximation. Iilhile the first approximation is certainly valid if the 
bands considered are broad, so that there exists a large energy separation between 
the valence bands and all the lower bands, the second approximation implies an 
averaging over angular variables, a procedure often used in the literature. 
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5. Donor Energy Level for Se in AS, Ga., As. 

In this section we like to apply some of the theory developed in the previous 

sections to a simple system. Selenium with the electronic structure [A] 3d^^4s^4p^, 

10 2 3 

when substituted for arsenic with the electronic structure [A] 3d 4s 4p in the 
crystalline compound AS, Ga As will evidently behave as a donor since it contri- 
butes one more 4p electron to the lattice. It is known experimentally that in 
pure GaAs (x=0), Se does indeed exhibit a donor level at 6 meV below the conduc- 


tion band edge. 

le c 
41) 


40) 


Recently, the position of this donor level within the band 


gap of the compound AS, Ga^ As as a function of x has been determined experi- 

X i-“X 


mentally.""^ With the aid of eq. (57) we shall now try to determine this donor 
level theoretically. In order to be able to do so, we must know the underlying 
band structure, E^(k) , and the interaction potential, Since we are dealing 

with donor states, the drastic simplification of ignoring all valence bands may 
be made. Furthermore, the conduction band of A2. Ga. As possesses one minimum 

X X*"X 

at the r point (k=0) and six equivalent minima at X ^(100), (100) ... (OOl)j . 

The relative position of those minima with respect to each other as well as their 
effective masses are a function of x, the mole fraction of A£,As in GaAs. Before 
we go into the details of the calculations, let us write eq. (57) again for the 
convenience of the reader, dropping the band index now that we are dealing only 
with the conduction band: 



i:)f(k) 


V(k - k')f(k") = 0 


k' 


Let us assume with Twose 


42) 


that we may write for f(k), 


f(k) = cj),(k - k.) , 
J ~ -J 


(71) 


( 72 ) 
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which is large only in the vicinity of a band extremum, such that 


E(k)f(k) 


E(kJ + (11^/2) 
“ j 


E-S 

a, 3 


.(k - k ) (k - k )„ 
3 j a j 3 




( 73 ) 


with sufficient accuracy. Here we have assumed a number of minima located at 
the positions k. and have performed a Taylor series expansion about k = k, retain- 
ing only the first two terms. Evidently 


3^E(k) 

8k^Bk3 




k=k. 


2m. - 

aa3 


(74) 


defines the effective mass tensor, in the usual way (a, 3 indicate Cartesian 

components). In the neighborhood of Ic = k^, eq. (71) now becomes 


E -1 


E(k.) - E + Cu^/2) 


a, 3 


(k - k.) (k - k.) 

a ^3 


. (k - k.) 


+ v(k - kO (j)^(k" - k^) = 0 

k'Ji 


(75) 


ilc * X* 

Multiplying eq. (75) by e' •- ~ and summing over all k yields 


e,k 




[e 


U 


E(k.) - (i,'/2)^ 


a, 3 


- kpe 


i(k-kj^) -r 


V(k - k')e^^— 
k" 


i(k" - l 5 ^) -r 


= 0 


(76) 
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By analogy to eqs. (55a) and (56), we define 


and 


F^(r) = 



i(k-kj*r 
^^(k - k^)e '' 


k 


’r 

V(k)e 

k. 


(77a) 


(77b) 


If we also define 


= E(k^) - + V(r) , (78) 

a, 3 

we may write eq. (76) in the C’jx<]pact form 

ik , • r 

e ^ (H, - E)F.(r) = 0 . (79) 

1 J ~ 

j 

The physical significance of the functions F^ may be seen if we retrace the steps 
which lead from eq. (32) to eq. (77a). Since is large only in the vicinity of 
k,, we may write 

F(k) = <j>,(k - k.) (80) 

Z— i 3 - ~3 

3 
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Then using eqs. (80), (77a) and (32) we find 






( 81 ) 


where are the pseudo Bloch functions associated with the conduction band 

minima located at and the wave function \l> of eq. (81) is an approximate 

solution of eq. (31). It follows now from both eqs. (79) and (81) that the 
energy can be expressed by 


E = 



(82) 


an expression which lends itself nicely to a variational calculation. But before 
we go ahead and do tills, we must delve Into the significance of eq . (79). 

Equation (79) constitutes one aquation for apparently many functions F^. Now, 
if we are dealing with a crystal like Si for instance, where only equivalent 
minima (6 in this case) exist, then all functions are equal by symmetry except 
for certain coefficients depending on the particular point group symmetry 
involved. But these factors are knov^m from group tlieory, and therefore eq . (79) 
is sufficient for a solution of the eigen-value problem. However in coses of 
non-equivalent minima, the argument fails and eq. (79) is not sufficient for find- 
ing a unique solution. Equation (79) has been derived from eq . (75) essentially 
by summing over alt k. But if the summation over k is restricted to a subzone 

centered about each minimum at k., we obtain a set of equations, one for each F,. 

'■J J 


36 


ORIGINAL PAGE 18 
OF POOR QUALITY 


This device has been used by Bassanl ot al- 


But this method leads also to 


non-local potentials in configuration space and gives significant simplifications 
only for potentials which are very strongly localized in momentum space. It 
would be therefore best to work with eq. (75), which is really a set of equa- 
tions, one for each rather than with eq. (79), which constitutes only 
one. equation for all Hov^ever we can avoid this complication if we use 

eq. (82) as a variational principle in wliich the function F^ belonging to non- 
equivalent minima are varied independently. Consequently, we set 


3 -r/a 

(ra^) e 


(83a) 


and for j = 1 to 6 


,,, , , 3. -r/b 

' i ~ ^ ^ 


(83b) 


representing the envelope functions belonging to the minimum at Y (83a) and the 

six equivalent minima at X (83b) of the conduction band of the A£^Ga^_^As. 

a and b constitute two independent variational parameters. The coefficients 

.t are assigned the value 1//6 for the ground state (the A., representation of 
J 

the group T^*"*^^). As model for the band structure, we take a simplified version 
of eq. (78). For the sole minimum at 1, we put 


«0 " - 4 ^ V(r) 
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and for the six equivalent minima at (k,0,0) (-k,0,0) (0,k,0) (0,-k,0) (0,0, k) 
and (0,0, -k) we. put 

2 

The zero of energy is measured from the position of the f minimum. For the 
potential V, wo use the semi-empirical expression (53). All quantities, the 
effective masses m^ and m^, the energy difference between direct and 

indirect band minima, k, the magnitude of the wave vector connecting the direct 
with any of the indirect minima, the strength of the attractive effective 
potential well and finally the static dielectric constant e(0) are functions 
of X. Values for the effective masses m^ and m^ and for the static dielectric 
constant t (0) have been used as reported by Hauser et al.^^\ Values for 
have been computed using band gap versus composition expressions culled from the 
same report. Continuing with the calculations, we insert eqs. (83) into 
eq. (82) and perform the indicated integrations. After some tedious but straight- 
forward algebra we obtain for E: 

E = + N 2 + + Ng + (85) 

where with 


c = ab/ (a + b) , 


D = 2 + 16/6 


|c/ (a + b)| 


3/2 


(1 + 


1 2 2, 
k c ) 


-2 


+ (1 + 


2 2 
kb) 


-2 


+ 16(2 + 


2 2 
Ic b ) 


-2 


(86a) 
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^2 

!3 - 

2 inya“ 


2 iUj^b 




C8bb) 


^2 “ ab 


r 

. _ 2 c . . / k . + .J, \ 

I + k^c'^ V" 0 ‘‘ "'I'V 


/,L . J- 

0 0 0 f 0 * 

(1 + 1^0*-)“ I rn^^a*^ nijb‘ 


,t?T,... 

2 injb' 


0 ') 

1 •}• 2 k^b" ^ 

"V + lb 


•) >1 2 


•H cf? 


(1 + kb) 


(1 + k^b^) + 16(2 + k“b“) -l- a/6 [c/(a-t-b)]^'‘-(l-t-k‘'c‘') 


‘} •) 
1 + k^'ir 


0 0 ' 
(2 -H iri)") 


-2 

•) 0 


2/1 


•) •) 


-•5 


(86c) 


= -c«/ 


1 - U -H 2 f + 2 \ ) 

i \ <• ' 

a 


(l + 2 « V,-2R/a 

ca \ a / 


(86tl) 


i 

0. 


N, = - >4 a(ab)“’^'^‘” ( \ c“(l + irc“) 


4 * 


2 /a:i 




+ lc“>-^) U - - lu) BlnltU + ixiskR 


/slnkR , ,„\| 

I + k c ' ' ) 


‘R/c 


-R/c 


N. - -.f/ I 


, + 0 k:\ ,- 2 r/'> 

I k - „ + - 


Bb 


1 4- 


R\ -2R/b 


(86f) 
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2, 2 /slnkR 


^6 “ ” ^ ^ \ kb 




2 . 2 


-2 


+ coskR - dvTCl + kb) 14- 2 co«2kR 


+ 


I - kb] t 


2 \kb 


kb) sln.2kR 


-2R/b^ 


, „./R ^ , 2, 2 /sinkR , ,,, ,,\ -2R/b 

+ (1 4- k b ) *1* t'OslkRj 0 , 


(Sbg) 


M ^ i2, 2 /sinv^kR . 

R? = - (2 + k b ) 4 - coa/2kRI e 

- A,r/(2 + k^b") ^1 - /2 ^K'OH/2kR + - kb^ alm'^2kb 

+ 8tt/ I (2 + k“b^) f/2 + oos/fkR'j 


-2R/b\ 


(8bb) 


2 

The laak now becomes formidable, (liven Che quantities R, i , k, e*" (the square 
of the electronic cliarge) and .*/, one must determine the variational paramelers a 
and b via the equations 


ha 


hb 


(B7) 


using of course, all expressions (86) in the defining equation (8h) for E. After 
having solved the transcendental eqs. (87) for a and b, the energy K final Iv is 
computed via eq. (8')). however, a and b signify the extent of the electron cloud 
for an electron bound to the impurity and, for shallow states, this extent is large, 
covering many lattice sites. The magnitude of k is of the order of a recipriuval 
lattice vector. Therefore it Is expected that the product ka or kh is large 
compai-ed to unity. If this is true, a drastic simplification arises since 
most of the terms of eqs. (86) become negligibly small. The neglect of the 
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terms containlnR ka or kb Is tantamount to the laoRU'ct of Intervaney mlxinp,. 


Nej.» lee ting then InUervalley mixing, the equation for the energy becomes 


^ + 




•f 


J 

^m^b* 


+ 1 
. , -r , j 


(Ng + 


( 88 ) 


where and N„ are given by eqs. (86d) and (86f), respectively, if we measure 

j ) 

all lengths in units of the Bohr I'adius a^j = 0.529A and energies in units of the 
Rydberg = 13.605/i eV, and if we introduce the quantities 

R = aajj, m^ ~ ^l"’e’ '’o “ 


where m is the mass of a free electron, the variational eqs. (87) become simply 
e 


y . 

V. e •’ = \ 

‘ .1 


2 

a (. 


'■'^/8jj)y; 


+ (2u/t ) (1 y i 



^ 0 or 1 . (90) 


The energy K becomes, with the aid of eqs. (90), 


i 

~n 




(91) 

Ror any given values of the parameters 'ij, a, i and eqs. (90) possess a unique 
set of solutions y . Once these solutions are found, the energy level K can be 
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computed via eq. (91). It is interesting to note that if c is made formally 
infinite, eqs. (go) become 


/j = yj (92) 

Letting t: >■ “ is tantamount to setting the Coulomb tail of the potential (53) equal 
to zero, and such a potential is representative of an isoelectronic impurity. But 
eqs. (92) only permit a solution (a bound state) if 

i e = 2.7183 , (93) 

J H 

and this inequality is almost identical to the well-known criterion for the 

onset of bound states derived from the one-band, one-site model for impurity energy 

45) 

levels by Koster and Slater. That theory leads to the expression 


1 = -^N 





(94) 


for the energy level of an isoelectronic impurity. If we now replace E^(k), 

2 2 St 

rather boldly, by the expression h 'k /2Ym and replace the 1 Brillouin zone by 

0 

a sphere of radius R = aa„, it can be shown that eq. (94) allows for a solution 
or energy level below the edge of the conduction band if the strength of the 
potential well satisfies the following inequality: 


Y« W/R^ ^ TT = 3.1416 


(95) 


However eqs. (90), valid for donor-like impurities with Z = 1 always possess 
a solution, or in other words, a bound state within the energy gap always exists. 
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Returning to the calculation of the energy elgen-value based on eqs. (90) and 
(91) , It Is necessary to obtain values for a and Remembering that a = 

and choosing for the cut-off radius R of the Coulomb potential (53) the 
covalent radius of Se, we find that a = 2.19. The last parameter to be deter- 
mined Is the strength of the attractive potential well ts^- ot eq. (53). A first- 
principle calculation via the defining eqs. (49) and (50) for the potential being 
a rather complicated undertaking, we choose to adopt a semi-empirical approach. 
Consequently, we put 




(96) 


and determined the values of the effective potential for pure A£As , 
and for pure GaAs, . , by means of the known energy values E % 6 meV in 
these two cases. In other words, taking the values for all relevant parameters 
like effective masses, etc., corresponding to either GaAs or A£As, we determined 
by means of eqs. (90) and (91) in such a manner that the energy E came out to 
have the experimentally determined value of 6 meV. In this manner we obtained 

-^GaAs - 'Si <57) 

Using these values and eq. (96) as well as the values of all the other param- 
eters as culled from ref. 44, we computed the energy eigen-value as a function of 

X. The results are plotted in fig. 1. Also plotted (as a dashed curve) is an 

41) 

average over experimentally determined values of the same quantity. Consider- 
ing the enormous simplicity of our theory, the overall agreement is rather good. 
The sharp maximum at the cross-over between direct and indirect band gap, experi- 
mentally 321 meV, is calculated as 316 meV. On the other hand the agreement is 
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rather poor in the "wings”, that is to say for either small values of 
X (x < 0.25) or for large values of x (x > 0.8), Here deviations of the com- 
puted values for the energy from the actual values are pronounced. One reason 
for this is due to the neglect of intervalley mixing. Another reason is due to 
the neglect of higher order terms in the expansion (eq. 73) of the energy E^C^)* 
This becomes clear when we look at the magnitude of the variational parameters 
a and b. It turns out that a varies betwen 207 (at x = 0) and 83 (at x = 1) 
Bohr radii and is thus comfortably large. This implies that the Fourier trans- 
form (eq. 72) is sufficiently peaked at k^ = 0 so that higher order terms in a 
Taylor series expansion about 0 may be neglected and eq. (73) is justified and 
at the same time, a mixing of this valley at k = 0 with the other six valleys at 
k. (j = 1 • • • 6) contribute negligibly to the energy. However b turned out to 
be much smaller than a. In fact, b varies between 2.68 (for x = 0) and 
2.35 (for X = 1) Bohr radii. In this case, the Fourier transform (72) becomes 
fairly delocalized, andeq. (73) constitutes a poor approximation. Also the 
product kb, where k is the magnitude of the distance (in k-space) between the 
central valley at the 1’ point and any of the six equivalent minima at X, now 
becomes small enough, so that intervalley mixing cannot be neglected. These 
shortcomings must be eliminated before an adequate theoretical understanding of 
the donor level of Se in A? Ga As becomes possible. 

X J.”“X 

For the sequel to this report, it is planned to investigate these topics 


further. 
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Fig. 1. Comparison of computed (dots) and experimental (dashed lines) 

values of the donor level of Se in Al^ ^^1-x ^ function 2 ^^ 

of X. The dashed lines arc averages over several measurements 
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